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Oh! Abstract 
D 

'■^ ' Mimicking pristine ^Z? graphene, we revisit the BBTW model for 4D lattice 

QCD given in [P. F. Bedaque et al. Phys. Rev. D78 (2008) 017502] by using the 
\ hidden SU (5) symmetry of the 4D hyperdiamond lattice Ti^. We first study the 

I link between the "^4 and SU (5); then we refine the BBTW 4D lattice action by 

^ I using the weight vectors Ai, A2, A3, A4, A5 of the 5-dimensional representation 

^ I of SU (5) satisfying ^ • Aj = 0. After that we study explicitly the solutions of 

I the zeros of the Dirac operator V in terms of the SU(5) simple roots ai, 0:2, 0:3, 

CX4 generating Ti/^; and its fundamental weights wi, uj2, ^3, ^4 which generate the 
reciprocal lattice H^- It is shown, amongst others, that these zeros live at the 
^ i sites of T-L^', and the continuous limit 2? is given by 7^k^ with d, 7^ and 

standing respectively for the lattice parameter of ?^4, the usual 4 Dirac matrices 
and the ^-D wave vector. Other features such as differences with BBTW model as 
well as the link between the Dirac operator following from our construction and 
the one suggested by Creutz using quaternions, are also given. 

Keywords: Graphene, Lattice QCD, 4D hyperdiamond, BBTW model, SU (5) 
Symmetry. 

1 Introduction 

In the last few years, there have been attempts to extend results on the relativistic elec- 
tron system on a 2D honeycomb (graphene) [HEKS] to a 4D honeycomb lattice (called 4D 
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hyperdiamond, denoted below as H4) and apply it to the lattice QCD simulations [1]-[T3]. 
These attempts try to construct Dirac fermion on 7/4 by keeping all desirable properties; 
in particular locality, chiral symmetry and the minimal number of fermion doublings 
[U E]; see also [12] and refs therein. In this regards, two remarkable approaches were 
given, first by Creutz suggesting an extension of graphene dispersion relations by using 
quaternions [U [8] ; and subsequently by Bedaque- Backoff- Tiburzi- WalkerLoud (BBTW) 
[5] proposing a 4D hyperdiamond lattice action with enough symmetries to exclude fine 
tuning. Apparently those two attempts look very close since both of them extend 2D 
graphene to 4D; however they have basic differences; some of them are discussed in [5]. 
The Creutz model involves a two parameter lattice action that lives on a distorted 4D 
lattice; and so looses the high discrete symmetry of the 4D hyperdiamond. The lattice 
action of BBTW model extends pristine 2D graphene; it is built on perfect 4D hyperdia- 
mond and has sufficient discrete symmetries for a good continuum limit. Nevertheless, in 
both Creutz and BBTW constructions, the distorted and perfect 4D hyperdiamonds are 
thought of as made by the superposition of two sublattices A4, and B4 with massless left 
and right-handed fermions as required by the no-go theorems for lattice chiral symmetry 

[laiis]. 

Guided by the rich symmetries of the 4D hyperdiamond ^4, we revisit in this paper, the 
BBTW model of ref [5] and its higher dimension extensions given in |T2] by using the 
hidden SU (5) [ resp. SU {d + 1)] symmetry of ^{4 [ resp. Tid+i ] and its reciprocal lattice 
Til [ resp. "H^+i ]. Focussing on 4D lattice QCD, we first review the link between BBTW 
construction and SU (5). Then we refine the hyperdiamond lattice action by using the 
weight vectors Ai, A2, A3, A4, A5 of the 5-dimensional (fundamental) representation of 
SU (5) as well as mimicking pristine 2D graphene which, in the language of groups, cor- 
responds precisely to SU (3). After that, we study explicitly the solutions of the zeros 
of the Dirac operator by using the SU (5) simple roots oci, ct2, ct^, ct^ generating H4; 
and its fundamental weights Ui, U2, W3, W4 generating the reciprocal lattice Til- We 
also comment the differences with BBTW construction; and exhibit the link between 
the Dirac operator, following from our approach, and the one suggested by Creutz using 
quaternions. 

The presentation is as follows: In section 2, we review briefly the BBTW parametrization 
of the real 4D hyperdiamond ^{4 and comment some particular discrete symmetries. In 
section 3, we study the link between Ha and the SU(5) symmetry. It is shown that 
T-L4 is precisely generated by the four simple roots cki, 0.2, cts, 0:4 of SU (5); and the 
reciprocal lattice HI is generated by its four weight vectors oji, 0^2, ^3, ^*^4,- In section 4, 
we revisit the BBTW model on 7^4 given in |5] and propose a refined 4D lattice action 
mimicking perfectly 2D graphene. In section 5, we study explicitly the zeros of the Dirac 
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operator; and in section 6 we re-derive the Borigi-Creutz fermions. In last section, we 
give a conclusion and make comments regarding other lattice models. 



2 4D hyperdiamond TLa 

Seen that the 4D hyperdiamond Ha, plays a central role in both BBTW and Creutz 
lattice models [H E], we start by studying this 4D lattice by exhibiting explicitly its 
crystallographic structure. In particular, we give the relative positions of the 5 first and 
the 20 second nearest neighbors and exhibit some particular discrete symmetries of . 
This analysis, which is useful in studying the link between the lattices and the SU (5) 
simple roots, is important in our construction; it will be used in section 3 to build the 
reciprocal lattice Til and in section 5 to study the dispersion energy relations as well as 
the zeros of the Dirac operator. 



In order to apply graphene simulation methods to lattice QCD, BBTW generalizes tight 
binding model of 2D graphene to the 4D diamond 7/4 O |6]; see also [3 El IS HO]. Like 
in the case of 2D honeycomb, this 4D lattice is defined by two superposed sublattices 

and B4 with the two following basic objects: 
First, sites in and B4 (L-nodes and R-nodes in the terminology of [5]) are parame- 
terized by the typical 4d- vectors rn with n = {ni,n2,n^,n^) and n^'s arbitrary integers. 
These lattice vectors are expanded as follows 



where ai, 3.2, as, a4 are primitive vectors generating these sublattices; and e is a shift 
vector as described in what follows. 

Second, the vector e is a global vector taking the same value V n; it is a shift vector 
giving the relative positions of the B4 sites with respect to the A4 ones; i.e: e = r'^ — rn, V 
n. In ref. j5], the a^'s and the e have been chosen as given by the following 4 - component 
vectors 



2.1 BBTW parametrization of Ti^ 



Aa 

S4 



Fn = ni ai + n2 ^2 + ^^3^3 + 714 a4 



r'n = Fn + e 



(2.1) 



ai = ei - 65 



as = 63 - 65 



e =65 



(2.2) 



a2 = 62 - 65 



a4 = 64 - eg 
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with the representation 



i (+75, +75, +75,+!) , 
i(+y5,-v^,-V5,+l) , 
i(-v^,-V5,+^,+l) , 



(2.3) 



i(-V5,+V5,-V5,+l) , 



and 



^1 ^2 ^3 "^4 



0, 0, 0, -1). 



(2.4) 



Notice also that the 5 vectors ei, e2, 63, 64, es define the first nearest neighbors to 
(0, 0, 0, 0) and satisfy the constraint relations. 



(2.5) 



COS 



4' 



« ^3 



showing that the e^'s are distributed in a symmetric way since all the angles "dij are equal 
to arccos (— |); see also figure (II]) for illustration. 

In the matrix representation (12. 3112. 4p . the free four vectors ei, 62, 63, 64 are permuted 
amongst each others by the typical unimodular matrices O^..^ acting as 



4 

u=l 



(2.6) 



[21] 



[32] 



(2.7) 



with, 

/ 1 \ 
0-100 

0-10 

V 1 / 

These transformations leave invariant the vector 
sub-symmetries of the permutation group generated by permutations of the five ej's. We 
also have 



/ -1 \ 

10 

0-10 

V 1 / 

= - (ei + e2 + 63 + 64); they are 



o 



o 



[43] 



[31] 



[32] 



o o 



o o 

[32]*^ [21] 

o o 

^[43]*-^ [31] 



(2.8) 



together with other similar relations. 
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2.2 Some specific properties 

From the figure ([T]) representing the first nearest neighbors in the 4D hyperdiamond 
and their analog in 2D graphene, we learn that each A4- type node at rn, with some 
attached wave function Aj.^, has the following closed neighbors: 5 first nearest neighbors 




Figure 1: On left the 5 first nearest neighbors in the pristine 4D hyperdiamond with 
the properties ||e,;|| = 1 and ei + 62 + 63 + 64 + 65 = 0. On right, the 3 first nearest in 
pristine 2D graphene with ||ej|| = 1 and ei + 62 + 63 = 0. 

belonging to with wave functions -Bm+dei ; and 20 second nearest neighbors belonging 
to the same A4, with the wave functions y4r^+d(e,-ej)- The first nearest neighbors are 
given by: 



lattice position 
Fn + dei 
Tn + de2 
Fn + de3 
Fn + de4 



attached wave 

Brn+de2 

Br^+dei 
-Brn+deg 



(2.9) 



Using this configuration, the typical tight binding hamiltonian describing the couplings 
between the first nearest neighbors reads as 



(2.10) 



where t is the hop energy and where d is the lattice parameter. Notice that in the case 
where the wave functions at Fn and Fn + de^ are rather given by two component Weyl 
spinors 

, B": 



^ n 




'dei 



Bi 



(2.11) 



ydei 
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together with their adjoints A'^^ and -Brn+de^' in the example of 4D lattice QCD to be 
described in section 4, the corresponding tight binding model would be, 



-t 



EE 



i=l 



+ he 



(2.12) 



where the e^'s are as in (12. 3 p and where the coefficients a'^^ will be specified later on. No- 
tice moreover that the term X^Li (^r„'^ad-^r„) vanishes identically due to Yl^=i — 0- 
The 20 second nearest neighbors read as 



rn±d (ei - 62) , rn±d (ei - 63) , rn±d (ei - 64) , 
rn±d (ei - 65) , rn±d (e2 - 63) , rn±d (e2 - 64) , 
rn±d (62 - 65) , rn±d (63 - 64) , rn±d (ea - 65) , 
rn±d (64 - eg) . 

At this order, the standard tight binding hamiltonian reads as follows 

5 



(2.13) 



rn i,i=l 

and in the case of Weyl spinors, we have 



(2.14) 



rn «,i=l 
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. ( A"- rr^ A°- -I- rr^ R"^ 

/ , \^^rn"ad^rn+a!(ei-e^) ^ ^r^" aa^rn+d(ei-ej) 



+ he 



(2.15) 



In what follows, we show that the 5 vectors ei, e2, 63, 64 65 are, up to a normalization 
factor namely precisely the weight vectors Aq, Ai, A2, A3, A4 of the 5-dimensional 
representation of SU (5); and the 20 vectors (ej — e^) are, up to a scale factor their 
roots = (Aj — Xj). We show as well that the particular property e^.e^ = which is 
constant V e^, V e^-, has a natural interpretation in terms of the Cartan matrix of SU (5). 



3 Link with SU (5) symmetry 

For later use, we exhibit here the hidden SU (5) symmetry of the 4D hyperdiamond; 



we show that 7^4 considered above is precisely the lattice Csu{b) studied in pJlj. More 
concretely, we show the three following: 

First, the 5 bond vectors ei, 62, 63, 64, 65 ( first nearest neighbors) are given by the 
5 weight vectors Ai, A2, A3, A4, A5 ( below, we set A5 = Aq) of the 5-dimensional 
(fundamental) representation of SU (5) which also satisfy 

Ao + Ai + A2 + A3 + A4 = (3.1) 
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We will show later that = with Aj.Aj = |. 

Second, the 4 primitive ones ai, a2, as, a4 used in generating Ti^ are particular linear 
combinations of the 4 simple roots cti, 0:2, 0:3, 0:4 of SU (5); see eq fl3.20p for the explicit 
relations. Recall that the SU (5) symmetry has 20 roots as given below, 

±0:1, ±(0:1 + 0:2), ±(01 + 02 + 03), ±(01 + 02 + 03 + 04) 

±02, ± (02 + 0:3) , ± (02 + 03 + 04) , 2^ 

±03, ± (03 + 04) , 

±03 

These vectors have all of them the same length o^ = 2; and so they generate the relative 
lattice positions of the second nearest neighbors in the 4D hyperdiamond. 
Third, the SU (5) has also discrete symmetries given by the so called Weyl group trans- 
formations generated by the cr^'s acting on generic roots (3 of SU(5) as follows, 

a«(/3) = /3-2^/3 = /3-(«./3)/3 . (3.3) 

These discrete transformations permute the roots (13.21) amongst themselves and are 
isomorphic to permutation group transformations. For instance, we have (oi) = 
-Oi and a^i (02) = oi + 02. 



3.1 Exhibiting the Hnk -Ha/SU (5) 

To exhibit explicitly the link between pristine lattice and the simple roots and the 
basic weight vectors of SU{5), we start by recalling some of its features; in particular 
the following useful ingredients: SU (5) is a 24 dimensional symmetry group; it has rank 
4; that is 4 simple roots Oi, 02, 03, 04; it has 20 roots ±/3jj given by eqf l3.2p . The 
simple roots Oi, 02, 03, 04 capture most of the algebraic properties of SU(5); and as a 
consequence, those of the 4D hyperdiamond crystal; in particular they generate the 20 
roots ±/3jj as shown on (13. 2 p and they have a symmetric intersection matrix K 
with inverse K^-^ given by. 



( 2 



-1 



V 












-1 



\ 



-1 
2 / 



K 





/4 


3 


2 


1\ 


1 


3 


6 


4 


2 


5 


2 


4 


6 


3 




V 1 


2 


3 


4/ 



(3.4) 



that encode the algebraic data of the underlying Lie algebra of the SU (5) symmetry. 
These simple roots define as well the 4 fundamental weights uj2, <^3, <^4 through the 
following duality relation 



(3.5) 
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These fundamental weights are important for us; first because they allow to build the 
reciprocal 4D hyperdiamond Til and second can be use used to expand any wave vector 
in Til as follows 

k = kiiVi + k2(jJ2 + ^3^3 + ^4^4- (3-6) 

From this expansion we read the relations fcj = k.CKj showing that the fcj's are precisely 
the wave vector components propagating along the CKj-directions; thanks to eqs fl3.5p . 



3.2 Other useful relations 

Using the matrices Kj^ and K^^-^, one can express the simple roots CKj in terms of the 
fundamental weight vectors LOf, and inversely the a;j's as linear combinations of the 
simple roots as given below, 



|q:2 + |q:3 + Ia4 



a;4 = icKi + |q:2 + |q;3 + 



(3.7) 



Using these relations, it is not difficult to check that they satisfy (13. 5p : for instance we 
have (jji.cti = | — | = 1 and u:i.ct2 = — | + | — | = 0; and similarly for the others u}2, 
LJ3, 0^4 and the intersections oji.ctj. Notice moreover that the fundamental weight vector 
oji defines a highest weight representation of SU (5) of dimension 5 with weight vectors 
Ao, Ai, A2, A3, A4 related to a; 1 as follows 



Ao 


= u:i 






Ai 


= U^i 


- «i 




A2 


= UJi 


- CKl 


— Ci2 


A3 


= 


- OCi 


- «2 - 0:3 


A4 


= iVi 


- OLl 


- ^2 - 0:3 - 


also 


express 


these vectors wei 




Ao 




^1 




Ai 




iV2-UJl , 




A2 




UJ3-UJ2 , 




A3 




^4-^3 , 




A4 




— a;4 



0:4 



(3.8) 



(3.9) 
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Furthermore, substituting uji by its expression (13. 7p . we get the following values of the 
Aj's in terms of the simple roots 



Ao = +|q;i + |q;2 + fas + |«4 
Ai = -|q;i + §0:2 + I0C3 + 
A2 = - + fcKs + 

A3 = -|q;i - |q;2 - fas + |a4 
A4 = -\ai - ^OL2 - |q;s - 



(3.10) 



5^ 5 

These weight vectors satisfy remarkable properties that will be used later on; in particular 
the three following: First, these Aj's obey the constraint relation X]i=o ~ ^ which 
agrees with f l3.ip and which should be compared with the identity 6^+62+63+64+64 = 0. 
Second, they have the intersection matrix 

Aj.Aj = I , Aj.Aj = -| , cos-(9ij = p^Tip^ = -J , (3-11) 

leading to eqf l2.5p . The third point concerns the zeros of the Dirac operator; see eqf l4.1ip 
to fix the ideas. They are given by solving the following constraint relations 



(3.12) 



where we have set 

Po = k.Ao , Pi = k.Ai , p2 = k.A2 
P3 = k.As , p4 = k.A4 

and where the phase ip = with an integer. The values of this phase are due 
to equiprobability in hops from a generic site at r to the 5 first nearest neighbors at 
r+^Y^Aj. This equiprobability requires 

5 

Y[e'^P^ = 1 = e''^ . (3.14) 
1=0 

Solutions of the constraint eqs fl3.12l) are then given by 

P^ = il7E ' ^ = 0,1,2,3,4 . (3.15) 

Notice moreover the two useful features: First, eqs( ]3.13l ) imply in turn that the wave 
vector k may be also written as 

k = poAo + PiAi + P2A2 + psAs + P4A4 (3.16) 

Multiplying both sides of this relation by Aj and using (13. lip , we find k.A^ = Pi — 
\ (po +P1+P2+P3+ Pi) = Pi] thanks to the identity Y.iPi = ^ following from ^ • Aj = 
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0. Second, expressing this vector k in terms of the basis u^i, 0^2, oj^, UJ4 of the reciprocal 
lattice; then using eqs fl3.9p giving the A/'s in terms of the Wi's, we get 



k = {po-pi)uJl + {Pl-P2)i^2 + iP2-P3)^^3 + {P3-P4)i^4 ■ (3-17) 

Putting back (13.151) . we find that the zeros of the Dirac operator are precisely located at 
the sites of the reciprocal lattice T-Ll- 

3.3 Link with BBTW parametrization of 

From eq fl3.8p . we can also determine the expression of the simple roots a^'s in terms of 
the weight vectors Aj's; we have: 

Oil = Ao — Ai , 0:3 = A2 — A3 , 
"2 = Ai — A2 , 0:4 = A3 — A4 

By comparing these equations with eq fl2.3ti2.3]) . we obtain the relation between the ej's 
used in [5] and the weight vectors of the fundamental representation of SU (5); 

Bi = -^Aj , Aj = ^^G?, 5 (3.19) 

Putting eqs (l3. 18113. 19p back into (12. 2p . we find that the 4 primitive vectors ai, a2, a3, a4 
generating the sublattice ^4 (resp. B4) are nothing but linear combinations of the four 
simple roots of SU (5), 



ai = -^oci 

a2 = (cKi + as) 

a3 = («! + 0:2 + 0:3) 

a4 = (CKi + 0:2 + 0:3 + 0:4) 



(3.20) 



From these relations, we read the identities of [5] 

These relations are just a property of Cartan matrix of SU (5). 
We end this section by giving the following summary: 

The 4D hyperdiamond ^4 is made of two superposed sublattices A4, and 84^. These 
sublattices are generated by the simple roots cki, cx.2, 0:3, cx.4 of SU{5). The relative 
shift vector between A4 and B4 is a weight vector of the 5-dimensional representation of 
SU{5). Each site in T-L^ has 5 first nearest neighbors forming a dimension 5 represen- 
tation of SU [5); and 20 second nearest ones; forming together with the "4 ^^ro roots", 
the adjoint representation of SU{5). The reciprocal space of the 4D hyperdiamond is 
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generated by the fundamental weight vectors uji, uj2, W3, UJ4 of SU (5). Generic wave 
vectors k in this lattice read as 



k = fciLJi + A;2<-J2 + k^ujs + fc4a;4 (3.22) 

where ki = {jpi-i — Pi) where pi is the momentum along the Aj-direction and (pi-i — Pi) 
the momentum along the CKj-direction in the real 4D hyperdiamond lattice "^4. In the 
particular case where all the momenta Pi = we have 

Y^Ke^''"^"' = e±*^(X^Af]=0 . (3.23) 
1=0 \i=o / 

This property will be used later on. 



4 BBTW lattice action revisited 
4.1 Correspondence 2D/4D 

To begin notice that a generic bond vector in Ha links two sites in the same unit cell 
of the hyperdiamond as shown on the typical coupling term ^r„-Bj^^^g.. This property is 
quite similar to the action of the usual 7^ matrices on 4D (Euclidean) space time spinors. 
Mimicking the tight binding model of 2D graphene, BBTW proposed in [5] an analogous 
model for 4D lattice QCD. There construction relies on the use of the following: First, 
the naive correspondence between the bond vectors and the 7* matrices 

Bi i — ^ 7i , i = l,...,5 , (4.1) 

with 

-es = 61 + 62 + 63 + 64 , 

-15 = 71 + 72 + 73 + 74 • (4-2) 



Recall that the four 7^^ matrices satisfy the Clifford algebra 7'^7'^ + 7^7^^ = 26^'^, 7^ = 
7^7^7^7^ gives the ± chiralities of the two possible Weyl spinors in 4D] and Fs is precisely 
the matrix F used in the Borigi-Creutz fermions [T7t [T8] : see also eq(2.5) of ref.[2T] for 
a rigorous derivation using SU (5) symmetry. Second, as in the case of 2D graphene, 
-4.4-type sites are occupied by left 0^. = (0r) right (pj^ = (^r^ 2-components Weyl 
spinors. i34-type sites are occupied by right Xr = (Xr+deJ and left Xl = (Xr+deJ Weyl 
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spinors. 





2D graphene 


4D hyperdiamond 


^4-sites at 






i34-sites at 
r„ + dBi 


R+ 


a 


couplings^ 


^rB^de. 
Br+de,A^ 


4 

4 



where the indices a = 1, 2 and a = 1, 2; and where summation over [i is in the Euchdean 
sense. For later use, it is interesting to notice the two following: In 2D graphene, the 
wave functions A^. and -Br+dei describe polarized electrons in first nearest sites of the 
2D honeycomb. As the spin up and spin down components of the electrons contribute 
equally, the effect of spin couplings in 2D graphene is ignored. In the J^D hyperdiamond, 
we have wave functions at each ^4-type site or i34-type one. These wave functions 
are given by: 

- the doublets 0"^ = (0r„,0r„) 0^ = (0r„'^r„) having respectively positive and 
negative 7^ chirality, these are the (^,0) and (0, |) representations of the SO ^ 
S\J (2) X SV (2). 

- the doublets Xr+de, = (xj+de,. Xr+deJ and Xr+de, = (xJ+de..X?+deJ havlug respectively 
negative and positive 7^ chirality. 

By mimicking the 2D graphene study, we expect therefore to have ^ kinds of polarized 
particles together with the corresponding "holes" as shown on the typical tight binding 
couplings 

(4.4) 

ef^li (xi+de, ^l) . ^i^22 (Xr+de, ^r) 

4.2 Building the action 

Following [5j, the BBTW action is a naive lattice QCD action preserving the symmetries 
of ?^4. To describe the spinor structures of the lattice fermions, one considers 4D space 
time Dirac spinors together with the following 7^ matrices realizations, 

7I = (g) , 72 = (g) or^ ^ ^3 ^ ^1 ^ ^3 ^ 

7^ = ® I2 , 75 = r3 ® I2 , ^ • ^ 

^ this correspondence differs from the one given by BBTW in [5]. 



12 




where the t"s are the Pauh matrices acting on the sublattice structure of the hyperdia- 
mond lattice I-L4, 

/ni\ / n ^ \ /in\ 

(4.6) 

The 2x2 matrices satisfy as well the Chfford algebra crV-' + a^a^ — 25'^^l2 and act 
through the coupling of left (f)^ (resp. 0^) and right Xr (resp. left Xl) 2-components 
Weyl spinors at neighboring A4- and B4- sites 



_(U _ _1 I a/5 _2 I V5„3 I IT 



where (7^ = (cr^, cr^, cr^, +il2) and = {a^ , , , —il2) ■ For later use, it is interesting 
to set 

_u _/i a/S _1 , a/S _2 , a/S _3 , i 

and similar relations for the other cr.ej and a.Si. 

Now extending the tight binding model of 2D graphene to the 4D hyperdiamond; and 
using the weight vectors Aj instead of ej, we can build a free fermion action on the lattice 
7^4 by attaching a two-component left-handed spinor 0" (r) and right-handed spinor 0^ 
to each ^4-node r, and a right-handed spinor ^ and left-handed spinor ^ 

to every i34-node at r + d^Xi. The action, describing hopping to first nearest-neighbor 
sites with equal probabilities in all five directions Aj, reads as follows:. 



4 



r i=0 

Clearly this action is invariant under the following discrete transformations 

'^%±d^x. <±.^A.(^j.): ' i^siT.K (4-10) 

Expanding the various spinorial fields in Fourier sums as J (^^^'^'^'Ck) with 

k standing for a generic wave vector in 1-L\, we can put the field action Sg^^^ into the 
form 



k 

where we have set 



4 4/4 



D ^Yl Ae*'"^^-^' = 5] I] Af e*'^^''-^' ) , (4.12) 

1=0 M=l V '=0 
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with 



Xj+tX^ Xf-zXf 



(4.13) 



and pi = k.Xi = Xluk/^Af. Similarly we have 



D =Y,Die 

1=0 



,1=1 



Mg-id^k.A, 



(4.14) 



j=o 



We end this subsection by making 3 remarks; the first one deals with the continuous 
limit; the second one regards the zeros of the Dirac operator and the third concerns 
the link with the Creutz fermions. In the continuous limit where the lattice parameter 
(i — )■ 0, we have 



E^ 

1=0 



^i.^±id^k.Xi 



E^r ± 



. 1=0 



+ 



(4.15) 



Moreover, since J2^=oK — ^^'^ because of the identity J2t=o-^f (k.A^) = k'^ following 
from eqsf l3.16ti3.17p . this limit reduces to 



E> 



1=0 



So we have 



D 



D 



(4.16) 



(4.17) 



fi=i 



The operators D and D have zeros for wave vectors k satisfying the following constraint 
relation 

k-A. = i0E > (4-18) 

with N an arbitrary integer. The point is that for these values, the phases e^'^^'^-^i = 
and the operators D and D get reduced to 



lj,=l \ 1=0 



lj.= l \i=o 



(4.19) 



which vanish identically due to the property J2t=oK ~ 0- Following [TT], the Dirac 
operator (14.111) in the Creutz lattice model reads as follows. 



z 

z* 



(4.20) 
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where z = 6qI + iOia^ + z^^20"^ + ^^30"^ with 

9i = sinpi + sinp2 — sinpg — sinp4 

62= sinpi - sinp2 - sinp3 + sinp4 (4 21) 

^3 = smpi — sinp2 + sinpa — sinp^ 

9q = B {AC — cospi — C0SP2 — cosps — C0SP4) 

and B and C two real parameters. In the Creutz lattice model, the zero energy states 
correspond to ^ = 0; this leads to the constraints 9i = which are solved by taking one 
of the momenta a.s pi = p and the others a.s pi = p or it — p. To make contact with our 
construction, the analogous of eqs fl4.2ip are given by: 



ei = 


Xle- 




+ X\e- 


id^Pi 


+ Xle- 


id^P2 


+ Xle- 


id^P3 


+ X\e 




02 = 


Xle- 


-id^po 


+ Xle- 


id^Pi 


+ Xle- 


id^P2 


+ Xle- 


id^PS 


+ Xle 




03 = 


Xle- 


id^po 


+ Xle- 


id^Pl 


+ Xle- 


id^P2 


+ Xle- 


id^PS 


+ Xle 


-id^Pi 


9o = 


Xte- 


id^fpo 


+ Xte- 


idyfpi 


+ Xte~ 


id^P2 


+ Xte- 


id^Pi 


+ X\e 


-id^P4 



(4.22) 



where pi = k.A;. These relations are complex and are, in some sense, more general than 
the Creutz ones fl4.2ip . The zeros of these solutions requires e = e**^ V Z = 0, 1, 2, 3, 4 
as anticipated in f l3.12p . 



5 Energy dispersion and zero modes 



To get the dispersion energy relations of the 4 waves components (f)^^, 0^, Xk 
their corresponding 4 holes, one has to solve the eigenvalues of the Dirac operator (14. lip . 
To that purpose, we first write the ^-dimensional wave equation as follows, 

J \ xk J V Xk 

where 0^ = (^ic'/'k)' Xk = (Xk' Xk) Weyl spinors and where the 2x2 matrices D, 
D are as in eqs (l4. 12114. 14p . Then determine the eigenstates and eigenvalues of the 2x2 
Dirac operator matrix by solving the following characteristic equation. 



(5.1) 



det 



( 





Dn 


D12 \ 





-E 


D21 


D22 




D21 


-E 





\ Di2 


D22 





-E J 



(5.2) 



from which one can learn the four dispersion energy eigenvalues Ei (k), E2 (k), E3 (k), 
£^4 (k) and therefore their zeros. 
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5.1 Computing the energy dispersion 



An interesting way to do these calculations is to act on (15. ip once more by the Dirac 
operator to bring it to the following diagonal form 



DD 
DD 



0k 



(5.3) 



Then solve separately the eigenvalues problem of the 2-dimensional equations DDip]^ 
E'^cj)]^ and DDx^ = E'^Xk- To do so, it is useful to set 



u (k) = ^^ + 1^^ , V (k) =^^ + id^ 



(5.4) 



with 



/i = l,2,3,4. 



(5.5) 



1=0 



Notice that in the continuous limit, we have 



u{k) 
vik) 



id^ (y + ik2) 
id^ (k3 + ik^) 



(5.6) 



Substituting (15. 4p back into (I4.12p and (14.140 . we obtain the following expressions, 



DD 



I |2 I I |2 
\u\ + |t>| 

2uv 



2uv 

\u\ + \v\ 



(5.7) 



and 



DD 



I |2 I I |2 

\u\ + \v\ 
2uv 



2uv 

I |2 I I |2 

\u\ + \v\ 



(5. 



By solving the characteristic equations of these 2x2 matrix operators, we get the 
following eigenstates (^^', Xk with their corresponding eigenvalues E^, 



eigenstates 



^k 



A 2/ 
^k 



2|«lbl 



^1 + 



2\u\\v\ 



2\u\\vl 



2\u\\vl 



eigenvalues 



El 



E'_ 



11 2 I 1 2 I I I I 

\u\ + \v\ +2 \u\ \v\ 



\uf + \vf — 2 \u\ If I 



(5.9) 



and 
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eigenstates 



2|u||t)| Ak \/ 2|M|k|Xk 



eigenvalues 



^2 ^ U|2 _^ |^|2 _^ 2 l^xl 



(5.10) 



^2 



\u\ + If I — 2 ImI If I 



By taking square roots of i?^, we obtain 2 positive and 2 negative dispersion energies; 
these are 

' (5.11) 



(5.12) 



E^ = +^{\u\±\v\r 

which correspond to particles; and 



El = -^{\u\±\v\ 
corresponding to the associated holes. 



5.2 Determining the zeros of E± and E. 



From the above energy dispersion relations, one sees that the zero modes are of two kinds 
as listed here below: 

zeros of both El = 0, = 
They are given by those wave vectors K.p solving the constraint relations u{K.f) = 
V (Kp) = which can be also put in the form 



for all values of /i = 1, 2, 3, 4; or equivalently like 



(5.13) 



(5.14) 



The solutions of these constraint equations have been studied in section 3; they are 
precisely given by eqs (13. 15113. 17p . Now, setting k = + q with small q = ||q|| and 
expanding D and D, eq (l5.ip gets reduced to 



af" 




Xk 



E 



Xk 



(5.15) 
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case ^2 = hut El = El^-^ ^ 
These minima are given by those wave vectors K = kmin solving the following constraint 
relation \u (K)| = \v (K)| or equivalently 



m,n=0 
m,n=0 



m,n=0 

with 



(5.16) 



Expanding this equahty, we get the following condition on the wave vector, 

4 

E Am COS (df K./3_) [tan (rff K./3_) - J^] = , (5.17) 



A - fX^X"^ _ Xl _ /'\3\4 _ xS n4\ 

•^nm V n m ^m^n) \^n^m ^m^n) 

» _ (\l \1 I \2 _ /x3 I x4 x4\ 

A possible solution is given by those wave vectors K obeying the relation K./3,, 
arctan (i3„m/-4-nm)- 



(5.1^ 



6 Re-deriving BC fermions 

In this section, we give the link between the above study based on SU(5) symmetry and 
the so called Borigi-Creutz {BC) model having two zero modes associated with the light 
quarks up and down of QCD. Recall that one of the important things in lattice QCD 
is the need to have a fermion action with a Dirac operator T> having two zero modes 
at points K and K' of the reciprocal space; so that they could be interpreted as the 
two light quarks. From this view, one may aslj^ whether there exists a link between the 
present analysis and the BC fermions [T71|T8]. In answering this question, we have found 
that the BC model can be indeed recovered from the analysis developed in this paper. 
In what follows, we give the main lines of the derivation. 



6.1 More on lattice action (14.91) 

One of the interesting lessons we have learnt from the analysis developed in the previous 
sections is that the lattice action for 4D hyperdiamond fermions may generally be written 

^we thank the referee for pointing out this question which allowed us to exhibit the relationship 
between our approach and BC fermions; see }21j for explicit details. 
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like, 

^ E ( E ^rT'^^r+aX, + '^r-aA. ) , (6.1) 

r \l=0 1=0 / 

where a = the weight vectors A/ as in eqsf l2.3ll3.19p and where and their complex 
adjoints are 4x4 complex matrices given by linear combinations of the Dirac matrices 
as follows 

E^'^m) ' (E^'^m) ' (6-2) 

,M=i / \m=i / 

with Q''^ linking the lattice euclidian space time index fi and the index I of the 5- 
dimensional representation of the SU (5) symmetry of the hyperdiamond. As such the 
lattice action (16. ip depends on the coefficients fij^ capturing 20 complex numbers that 
form a 5 X 4 matrix representing the bi-fundamental of SO (4) x SU (5) 



/ fio Q\ Ql Ql Qf \ 

(6.3) 



^"^2 ^^2 ^^2 
^"^3 ^^3 ^^3 

V Ql Ql Ql Ql Qj j 



This rank two tensor, which we decompose as (w^, fi^) with = fij^ a complex 4 
component vector and fij;^ a complex 4x4 matrix, gives enough freedom to engineer 
Dirac operators with a definite number of zero modes. Below, we derive the constraint 
equations for the zero modes of the Dirac operator; and in next subsection we apply the 
analysis to the BC model. 

6.1.1 Dirac operator 

In the reciprocal space, the lattice action (16. ip reads as 

5~ EiE'^kl^^k) (6.4) 

k \At=l / 

with Dirac operator reading as follows 

^ = iE^'(^M + ^.), (6.5) 
/i=i 

and where and its complex adjoint -D^ are given by: 

4 4 

D^, = Y^\e^''^■^^ , 5^ = ^ nj^e-^""^-^' . (6.6) 

1=0 1=0 
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These operators depend on 40 = 2 (4 + 16) real numbers 

uj, = \{u, + iv,) , n-;^ = iR- + ij;: , (6.7) 

and also on the five momenta pi = hki along the A/- directions. Since ki = k.A/ and 
because of SU(5) symmetry we have moreover the constraint relation 



ko + ki + k2 + k3 + k4 = 0, mod ^ , (6.8 



allowing to express one of the five ki's in terms of the four others. For instance, we can 
express k^ as follows: 

2n 

ko = - (ki + ^2 + ^3 + h) , mod — . (6.9) 

a 

The next step is to find the set of the wave vectors k^ = {ki, k2, k^, k^) that give the zeros 
of the Dirac operator. These zeros depend on the numbers m^, v^, i?^ and which can 
be tuned in order to get the desired number of zeros. 

6.1.2 Zeros modes 

The zero modes of the Dirac operator V given by eqs (16.5116^ are obtained by solving 
the following constraint equations 

4 4 4 4 

(^M + cos ak + tJ2J2^' K - sin ak = 0, (6. 10) 

At=l 1=0 /i=l 1=0 

together with the constraint eq fl6.8l) . Using the decomposition Q'-^ = (a;^,^^), we can 
decompose these constraints as follows 

A / A 4 \ 

fj.= l \u = l u=l J 

where we have set 

4 4 

A = cos akQ (w^ + Q^) + i sin ak^ [u^ — Q^) . (6.12) 

Moreover using (16. 7p we can put the above constraint relations into the following equiv- 
alent form 

4/4 4 \ 

fj,= l \u=l u = l J 

and 

A = cosafco I'^u^ - sinafco 7^^^^ = 0, (6.14) 
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with ko given by eq fl6.8p . Eqsf l6.13ti6.14p define a highly non hnear system of coupled 
equations in the four fcj,'s; and are difficult to solve in the generic case. To overcome 
this difficulty, one may deal with these equations by focusing on adequate solutions for 
the kjj^s; and engineer the corresponding tensor. Below, we apply this idea to the BC 
model. 



6.2 BC fermions 
6.2.1 Deriving the model 

Borigi-Creutz model [17] is a simple lattice QCD fermions for modeling and simulating 
the interacting dynamics of the two light quarks up and down. The Dirac operator of 
this model reads in the reciprocal space as follows, 

4 4 4 

'^BC ~ ^^7^sinafc^-^^7''cosafc^ + ^^rcosafc^-fr , (6.15) 

11=1 IX=1 /i=l 

with r = I (7^ + 7^ + 7^ + 7^). From this expression, one can check that this operator 
has two zero modes given by the two following wave vectors, 

(1) : ik,,k2,ks,k,) =(0,0,0,0) , ^g^g^ 

(2) : {ki,k2,k3,ki) =(^,^,^,^) , 

satisfying the remarkable property 

27r 271 

ki + k2 + ks + k4 = — , mod — . (6-17) 
a a 



Clearly the operator Vbc corresponds to a particular configuration of the complex tensor 
Q"^^ and the vector u^. To see that is indeed the case, notice first that the matrix F can 
be conveniently rewritten as F = |'i9^7'^ with, 

= (1,1,1,1) • (6.18) 

The same feature is valid for the sum Y2^=i cos cik^, which can be also put in the form 
St=i cos afcj,. Putting these expressions back into the above Vbc relation, we get: 

'Dbc ~ - V YK^^^aK - - V 7''M;;cosafc^ - -F, (6.19) 

with 



a ^ — ^ a ^ — ^ a 

fM,l'=l fl,l/=l 
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or more explicitly, 



/ I 1 _i _i _i 

' ^2 2 2 2 

1 ,1 _1 _1 

'2 ~^ 2 2 2 

1 _1 ,1 _1 
'2 2 ~'~2 2 

2 2 2 ^2 / 



V 



(6.21^ 



Now comparing eq( l6.19p with the general form of the Dirac operator of eqf l6.5ll6l6|) which 
also reads like 



4/4 4 

A + ^ i^R"^ cos aK - l^Jl sin aK ) = 0, 

^=1 \u=l 



(6.22) 



we see that Vbc can be recovered by taking, 



(6.23) 



and 

A = -2r = - (7I + 7^ + 73 + 7^) . (6.24) 

Eq (l6.23p leads to il'^ = (M^ + ; by substituting by its expression given above, 
the tensor Q'j^ reads more explicitly like 



(6.25) 



The second constraint relation f l6.24p requires 

{u^ cos afco — f sin ako) 



(6.26) 



with as in fl6.18p . Moreover, using eqs fl6.17f6.9p . we end with 



(6.27) 



and a free vector which, for simplicity, we set to zero. Thus the tensor fi^^ = (cj^, fij^) 
describing the BC fermions is given by 



/-I 



l+2i 1 



v 



4 



1 1 
'2 4 



1 1 
'2 4 



1 

4 4 
l+2t 1 
4 4 



l+2)i 1 
4 4 
1 l+2t 



(6.28) 



with the trace property = "I'^^m- 
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6.2.2 Symmetries 

Here we want to make a comment on particular symmetries of 4D lattice QCD fermions 
by following the analysis of ref . |6] where the study of the renormalization of this class of 
models has been explicitly done. There, it has been found that the breaking of discrete 
symmetries, such as parity V : \E'(A;, k^) — )■ 'y^'^{—k, /C4) and time- reversal T: {k, k^) — )• 
7574\E'(fc, — /C4), is behind the appearance of relevant dimension 3 operators O^^^ and 
marginal dimension 4 ones O^^ in the analysis of the Symanzik effective theory with 
lagrangian £e// = ^„ a" c£^On^ . Following the above mentioned work, one starts 
from the 4D lattice action, 

1 4 

^x+a., - K A' - 2iBC V1/+ 7^ V]/ J . (6.29) 

which depends on two real parameters B and C that are fixed by physical requirements 
and symmetries. This typical action depends also on particular combinations of gamma 
matrices A'^ = X]^=i l'^ where the coefficients A'^^ given in [6], form an invertible 4x4 
matrix with det (^t') — —16^5. Notice that setting S = 1, C = ^ one recovers the 
Borigi action. By performing transformations of fl6.29p using Fourier integrals to move 
to the reciprocal space, similarity operations to exhibit particular symmetries, expansion 
in powers of the lattice spacing parameter a to use the Symanzik effective theory; and 
switching on the usual gauge interactions — )■ = — igA^ with field strength 
= ^ l^fi, '^u], we get up to the first order in the parameter a the following effective 
field action, 



Q ® /) D^Q - ^ j;.. J-^" + aCs + ord (a^) 



with Qa standing for the quark isodoublet {ua,da) and O5 some dimension 5 operator 
that can be found in |6]. This effective theory has several symmetries in particular: 
(1) the manifest gauge invariance, (2) Ub (1) baryon number, (3) Ul (1) x Ur (1) chiral 
symmetry, (4) CVT invariance and (5) symmetry under 5*4 permutation of the four 
hyperplane axis corresponding to C = BS with e*^ = C + iS. The authors of [6] 
concluded their work by two remarkable results: (i) the engineering of chirally symmetric 
action with minimal fermion doubling which does not generate dimension three operators 
is possible as far as VT symmetry is preserved. This invariance is sufficient to forbid 
the relevant dimension 3 operators O^^ whose typical forms are listed below, 

broken P : O^'^ = ^ 7, ^^fc^, , (^f = ^ ^475 

broken r : = 2 7,75 , = ^ 7,75 
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with 7j standing for 7^^, 72 , 73. (ii) For particular values of parameters of the theory, 
there may emerge some additional non standard symmetries which could be used to 
eliminate the relevant operators. These results are important and may serve as guide 
lines in dealing with this problem by using the hyperdiamond symmetries based on roots 
and weights of SU (5). Below, we give a comment on this matter; an exact answer needs 
however a deeper analysis. In the SU (5) framework, the previous action (16.291) gets 
extended as follows 

'^-^ ~ ^ E E ( E [^^l^^^M, + ^.Y^l^.-ax,] J , (6.30) 

X \ i=0 / 

where as before and the A^'s are the weight vectors of the 5-dimensional representation 
of SU(5). Clearly, this lattice action is more general than eqf l6.29p : and has two inter- 
esting features that are useful in dealing with the study of underlying symmetries and 
renormalization of Sgu^- First, the SU {5) property fl3.ip on the weight vectors namely 
Xf = induces in turns the following constraint relation on the wave vectors A;^, 

5 5 
^ fc; = 0, with ki = ^ k^.Xf. 

1=0 ^i=l 

This constraint is invariant under VT symmetry acting on wave vectors as — )■ —k^; but 
not preserved under parity V nor time- reversal T separately. Second the generalized 
action Sg^^ depends on 20 complex {40 real) moduli carried by the tensor Q'-^. This 
number gives quite enough freedom to engineer QCD-like models with two zeros for the 
Dirac operator as we have done in case of BC model; may lead to desired symmetries 
of the Symanzik effective theory that follow from the expansion of the action Sgur, in 
powers of the lattice parameter; and may allow to make appropriate choices to eliminate 
relevant operators. Progress in this matter will be reported in a future occasion. 

7 Conclusion 

In this paper, we have studied the lattice fermion action for pristine 4D hyperdiamond 
7/4 with desired properties for 4D lattice QCD simulations. Using the SU (5) hidden 
symmetry of ^{4, we have constructed a BBTW- like lattice model by mimicking 2D 
graphene model. To that purpose, we first studied the link between the construction of 
[5] and SU (5); then we refined BBTW lattice action by using the weight vectors Aq, Ai, 
A2, A3, A4 of the 5-dimensional representation of SU (5). After that we studied exphcitly 
the solutions of the zeros of the Dirac operator in terms of the SU (5) simple roots cki, 
CX.2, CKs, 0:4, and its fundamental weights oji, LO2, ^*J3, ^^i- We have found that the zeros 
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of the Dirac operator live at the sites k =-^^ (^i^i + N2LO2 + N^oj^ + N4UJ4) of the re- 
ciprocal lattice T-tl; with Ni integers. In addition to their quite similar continuum limit, 
we have also studied the link between the Dirac operator following from our construction 
and the one suggested by Creutz using quaternions; the Dirac operator in our approach 
may be viewed as a " complexification" of the Creutz one where the role played by the 
sin(pj)'s and the cos(pj)'s is now played by e*^' as shown in eqs fl4.2ip and fl4.22p . The 
exact link between our approach and the Borigi-Creutz fermions has been worked out 
with details in section 6; where it is shown that the BC action follows exactly from (16. ip 
with eqs fl6.2ll6.28p giving the linear combinations of the Dirac matrices of the model. 




It is also interesting to notice that our approach is general; and applies straightforwardly 
to lattice systems in diverse dimensions. The fact that the 4D hyperdiamond is related 
to SU (5) fundamental weights Wi, u}2, ^3, ^4, and its simple roots cti, cx.2, ct^, 0:4 is not 
specific for 4-dimensions; it can be extended to generic dimensions D where the underly- 
ing D-dimensional hyperdiamond lattice has a hidden SU {D + 1) symmetry with simple 
roots Ck.1,..., OL£) and fundamental weights cji, uJd- From this view, the 2D graphene 
has therefore a hidden Sf/ (3) symmetry as reported in details in [3J. Our construction 
applies as well to the fermion actions given in [T2] . 
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